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We report measurements of the momentum dependence of the superconducting gap in Bi^Sr^C 'aCu20s+x 
(Bi2212) with angle-resolved photoemission spectroscopy using a dense sampling of the Brillouin zone in 
the vicinity of the Fermi surface. In the Y quadrant of the zone, where there are no complications from 
ghost bands caused by the superlattice, we find a gap function consistent within error bars to the form 
cos(k x ) — cos(ky) expected for a d-wave order parameter. Similar results are found in the X quadrant with 
the photon polarization chosen to enhance main band emission over that due to ghost bands. 

PACS numbers: 74.25.-q, 74.72.Hs, 79.60.Bm 
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The past few years have seen an intense effort on the 
part of the condensed matter physics community to de- 
termine the form of the order parameter in the high tem- 
perature cuprate superconductors!! Part of this effort 
has relied on the use of high resolution angle-resolved 
photoemission spectroscopy (ARPES) since this is the 
only known technique which can directly map out the 
momentum dependence of the superconducting gapBtJ. 

Very recently there has been considerable progress 
in ARPES studies of the cuprates. This includes evi- 
dence for a spectral function interpretation of ARPES 
via the^tudy of sum rulesa, a detailed, study of the Fermi 
surfaced and momentum distributiorfl in Bi2212_and an 
elucidation of the polarization selection rulesEHa in the 
presence of the structural superlattice, which necessar- 
ily affects ARPES data on the highest quality Bi2212 
samples. All these points are crucial in accurately deter- 
mining the momentum dependence of the gap. 

Modeling the experimental data in terms of soectral 
functions gives a quantitative estimate of the gapa at a 
given location on the Fermi surface, but perhaps even 
more importantly, it provides a framework within which 
one can think about the data. This points to the im- 
portance of having a very fine momentum space grid for 
data, and also for a very precise determination of the 
Fermi surface, in the neighborhood of the region .where 
nodes in the gap occur. Further, the complicationsa aris- 
ing from the superlattice make it—simplest to study the 
gap anisotropy in the Y quadrantL9, since here the main 
CuO Fermi surface is widely separated fcom the "ghost" 
Fermi surfaces due to the superstructure!! 

In this paper, we summarize the results obtained on 
several samples of Bi2212 with experiments which were 
designed taking into account all of the points listed 
above. We find that, within error bars, the excitation 
gap on the main Fermi surface sheet is consistent with 
cos(k x ) — cos(fc y )|. The ARPES experiment cannot of 



course measure the phase of the order parameter, but this 
result is strongly suggestive of a d x 2_ y 2 order parameter. 
Finally, in order to prove that the two node gap previ- 
ously seen in the X quadrant!] was indeed an artifactEl of 
the superstructure, we also study the X quadrant using 
a polarization geometry which enhances the main band 
relative to the "ghost" band, and obtain results which 
are in agreement with the Y quadrant. 

The experiments were done at the Synchrotron Ra- 
diation Center, Wisconsin, using a high resolution 4-m 
normal incidence monochromator with a resolving power 
of 10 4 at 10 11 photons/s. We used 22 eV photons, with 
a 17 meV (FWHM) energy resolution, and a momentum 
window of radius 0.0457T (in unitst^l of 1/a*). 

The crystals, which were grown by the traveling sol- 
vent floating zone method with an infrared mirror fur- 
nace, have low defect densities with structural coherence 
lengths ~ 1250 A obtained from X-ray diffraction. The 
samples were cleaved in-situ at 13 K in a vacuum of < 
5xl0 -11 Torr. Most samples have very flat surfaces after 
cleaving, as measured by specular laser reflections. A flat 
surface is crucial for this experiment since it directly af- 
fects the momentum resolution. Another measure of the 
sample quality is the observation of the "ghost" bands 
due to the superlattice distortion; in our best sample 
(87K T c , with a IK transition width) we have now also 
seen evidence for the second umklapps from the superlat- 
tice. We will also present data from three other samples 
with a 90K T c . 

In Fig. 1, we show the 13K experimental energy distri- 
bution curves (EDCs) for the 87K T c sample for various 
locations on the main band Fermi surface (FS) in the Y 
quadrant. The spectra shown correspond to the mini- 
mum observable gap along a set of k points normal to 
the FS (for a detailed discussion of this in the context of 
a study of particle-hole mixing, see Ref. |J). These spec- 
tra are obtained from a dense sampling of k-space in the 
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vicinity of the FS which is almost five times denser than 
previous data. In Ref. [5] a k-step size of 0.0647T normal 
to the FS and 0.064-0.081tt along the FS was used in the 
Y quadrant (the radius of the k- window is 0.0457r). The 
problem is that the bands are highly dispersive along TY 
with an energy change of about 85 meV per 0.064-7T step, 
and thus k^ cannot be located accurately. Therefore, 
in our new measurements, we use a step size of 0.02257T 
normal to the Fermi surface and 0.0457T along the Fermi 
surface. This not only allows us to map out the nodal re- 
gion more clearly, but also improves by a factor of about 
three our ability to locate that spectrum whose binding 
energy at the center of the momentum window is closest 
to the Fermi energy (that is, the step size normal to the 
Fermi surface now corresponds to a dispersion of about 
30 meV per step). In addition, at each k point the pho- 
ton polarization is chosen along TX so as to maximize 
emission on the TY diagonal direction, i.e., TY _L geom- 
etry (a polarization rotated 45° relative to this was used 
in Ref. g). 

In each panel of Fig. 1 we also plot the spectrum of 
a platinum reference, in electrical contact with Bi2212, 
measured periodically to determine the chemical poten- 
tial, and to check for drifts in beam energy. (The poly- 
crystalline Pt spectrum is a weighted density of states 
whose leading edge is an energy-resolution limited Fermi 
function). The simplest gap estimate is obtained from 
the mid-point shift of the leading edge of Bi2212 rela- 
tive to Pt. This has no quantitative validity (since the 
Bi2212 EDC is a spectral functiorfj, while Pt is a den- 
sity of states) but yields an angular dependence which is 
qualitatively similar to the results described below. 

The simplest way to quantitatively estimate the gap is 
to modelD the data in terms of a simple BCS lineshape 
formula, taking into account the measured energy disper- 
sion and the known energy and momentum resolutions. 
Two important points need to be discussed in connec- 
tion with such fits: first, the lack of knowledge about 
the spectral lineshape, especially its incoherent part, and, 
second, the large background in Bi2212 whose origin is 
unclear. In the large gap region near the M point, we 
see a linewidth (imaginary self-energy) collapse, for fre- 
quencies smaller than 3A, upon cooling well below Tjj'Lj. 
Thus the coherent piece of the spectral function is mod- 
eled by the BCS lineshape, with all of the incoherent 
part lumped together with the experimental background 
in so far as the fitting procedure is concerned. We also 
showedO that it was self-consistent (though perhaps not 
unique) to make the same assumption in the small gap 
region; the much larger width of the EDC in the diagonal 
direction arising due to the k resolution combined with 
the large dispersion. 

The gaps extracted by fitting the spectra of Fig. 1 are 
shown in Fig. 2. We emphasize that since the gap is de- 
termined by fitting the resolution-limited leading edge of 
the EDC, it is relatively unaffected both by self-energy 
effects, and by the experimental background which cuts 
off at low frequencies. To check this, we have made an 



independent set of fits to the small gap data where we 
do not use any background fitting function, and only try 
to match the leading edges, not the full spectrum. The 
two gap estimates are consistent within a meV. The (ver- 
tical) error bars of ± 3meV on the gaps in Fig. 2 come 
primarily from the quality of the fit to the leading edge of 
the data and uncertainty in the location of the center of 
the k- window (near the diagonal direction) , with smaller 
contributions coming from chemical potential determina- 
tion (± 0.5meV) and background modeling (± 0.5meV). 
Horizontal error bars represent the accuracy in which we 
can determine the Fermi surface angle. This does not 
include any effective error bar coming from k resolution, 
since this is in principle taken into account in the fits. 
The angular variation of the gap in Fig. 2 is in excel- 
lent agreement with that expected from_a d-wave order 
parameter of the cos^) — cos(k y ) forn£3. 

Next we turn to the X xjuadrant of the 87K T c sam- 
ple. It is now recognizee!] that the two node gap ob- 
served previously^ in the TX\\ geometry came from an 
even linear combination of the "ghost" bands due to the 
superlatticeO. We thus choose the polarization along 
TY, so that on the diagonal we are in a TX J—geometry 
which enhances emission from the main bandlij. The X 
quadrant gaps, determined from spectral function fits, 
are plotted in Fig. 3. We see thai, the hump in the 
gap along TX (45°) seen previouslyB has indeed disap- 
peared. The solid curve is a fit of the data to a d-wave 
gap function with a small sample misalignment (1.4° in 
real space). Note that for this data set, the step size 
along the Fermi surface was 0.1357T and so is not dense 
enough around 45° to address the question of the detailed 
behavior around the node. 

We next summarize the results (Fig. 4) for Y quad- 
rant gaps extracted from fits on three different 90K T c 
samples. The main point is to note possible complica- 
tions which arise in interpreting data sets which are not 
as dense in k space as the detailed Y quadrant mea- 
surements on the 87K sample described above (the step 
size along the Fermi surface being twice as large). The 
results on sample I have a region of reduced gap, consis- 
tent with zero, near 45°. To some extent this may be an 
artifact of thcpiinitc diameter of the k-window, which is 
6° in FS angleEil. In addition, we found that for the small 
gap points of sample I the leading edge of the data lies 
above, i.e., to the right of, that of a zero gap spectral fit, 
assuming the k-window was centered at kp. We find that 
a combination of factors (k-window center, chemical po- 
tential drift, and background) discussed above can indeed 
conspire to produce such an anomalous shift. These fac- 
tors are already reflected in the ±3 meV error bars on the 
gap estimate, but we reiterate that the error bars must 
be borne in mind while interpreting the final results. 

Similar issues arise with the results on sample II, and 
in addition both this and sample III appear to have an 
asymmetry about the diagonal (45°). For sample III, 
the asymmetry is quite large and almost certainly due to 
sample misalignment (of 1.8° in real space); the sample 
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Ill results in Fig. 4 have been shifted by the correspond- 
ing FS angle to facilitate comparison with the other two 
samples. All of these caveats would also apply to the 
sparse Y quadrant data presented in Ref. [|, and, in par- 
ticular, the shift of the node away from 45° in that data 
may also have been a consequence of sample misalign- 
ment. 

In summary, we would say that the results of the 90K 
samples are consistent, within error bars, with those ob- 
tained from the detailed measurements on the 87K sam- 
ple. There is no hard evidence for cither an extended 
region of nodes about 45°, or for a mixed symmetry gap. 

In conclusion, we have determined the intrinsic mo- 
mentum dependence of the superconducting gap in 
Bi2212 using dense momentum sampling of the zone in 
the vicinity of the Fermi surface. We find that the olt. 
served gap is consistent with the simple d-wave formt£l 
cos(fca;) — cos(fc y ) within experimental error bars. 
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as the main band, should also epntribute in this geometryu. 
There is, however, no evidenceQ in the data for such an odd 
"ghost" band, possibly due to final state effects. We note 
that the emission intensities from the main and ghost bands 
are not well understood at this time. 

14 So as not to prejudice the fits, the gap was considered to be 
a constant within the momentum window. This condition is 
obviously violated in the d-wave case near the node where 
the gap is varying rapidly with momentum. 

15 Such an order parameter arises naturally from theories with 
antiferromagnetic spin fluctuations and/or strong correla- 
tions. For recent reviews, see D. Pines and P. Monthoux, 
J. Phys. Chem. Solids 56, 1651 (1995) and D. J. Scalapino, 
Phys. Rep. 250, 329 (1995) and J. Phys. Chem. Solids 56, 
1669 (1995). 

FIG. 1. Bi2212 spectra (solid lines) for a 87K T c sample 
at 13K and Pt spectra (dashed lines) versus binding energy 
(meV) along the Fermi surface in the Y quadrant, with loca- 
tions shown in Fig. 2. 

FIG. 2. Y quadrant gap in meV versus angle on the Fermi 
surface (filled circles) with fits to the data using a d-wave gap 
(solid curve). Labels of data points correspond to the spectra 
of Fig. 1. Inset shows their locations in the zone as well as 
the photon polarization direction. 
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FIG. 3. X quadrant gap in meV, for the 87K T c sample, 
measured at 13K, versus angle on the Fermi surface (filled 
circles) with fits to the data using a d-wave gap (solid curve) . 
The photon polarization is along FY. 



FIG. 4. Y quadrant gap in meV, measured at 13K, versus 
angle on the Fermi surface for three different Bi2212 samples 
each with a 90K T c . For visual clarity only a representative 
error bar has been shown. 
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